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Abstract 



We study the resummation of large logarithmic perturbative corrections to the single- 
inclusive jet cross section at hadron colliders. The corrections we address arise near the 
threshold for the partonic reaction, when the incoming partons have just enough energy 
to produce the high-transverse-momentum final state. The structure of the resulting 
logarithmic corrections is known to depend crucially on the treatment of the invariant 
mass of the produced jet at threshold. We allow the jet to have a non-vanishing mass 
at threshold, which most closely corresponds to the situation in experiment. Matching 
our results to available semi-analytical next-to-leading-order calculations, we derive 
resummed results valid to next-to-leading logarithmic accuracy. We present numerical 
results for the resummation effects at Tevatron and RHIC energies. 
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1 Introduction 



High-transverse-momentum jet production in hadronic collisions, H1H2 — ► jet X, plays a fun- 
damental role in High-Energy Physics. It offers possibilities to explore QCD, for example the 
structure of the interacting hadrons or the emergence of hadronic final states, but is also inti- 
mately involved in many signals (and their backgrounds) for New Physics. At the heart of all 
these applications of jet production is our ability to perform reliable and precise perturbative cal- 
culations of the partonic short-distance interactions that generate the high-transverse-momentum 
final states. Up to corrections suppressed by inverse powers of the jet's transverse momentum pt, 
the hadronic jet cross section factorizes into parton distribution functions that contain primar- 
ily long-distance information, and these short-distance cross sections. In the present paper, we 
address large logarithmic perturbative corrections to the latter. 

At partonic threshold, when the initial partons have just enough energy to produce the high- 
er jet and an unobserved recoiling partonic final state, the phase space available for gluon 
bremsstrahlung vanishes, so that only soft and collinear emission is allowed, resulting in large 
logarithmic corrections to the partonic cross section. To be more specific, if we consider the cross 
section as a function of the jet transverse momentum, integrated over all jet rapidities, the par- 
tonic threshold is reached when y/s = 2p?, where \/s is the partonic center-of-mass (cm.) energy. 
Defining xt = 2p T j y/s, the leading large contributions near threshold arise as a^ipr) ln 2m (1 — x^) 
at the /cth order in perturbation theory, where m < k (the logarithms with m = k are leading) 
and as is the strong coupling. Even if px is large so that «s(pt) is small, sufficiently close to 
threshold the logarithmic terms will spoil the perturbative expansion to any fixed order. Thresh- 
old resummation [H [21 [31 HJ El EJ [7] , however, allows to reinstate a useful perturbative series by 
systematically taking into account the terms a§ ln 2m (1 — x T ) to all orders in as- This is achieved 
after taking a Mellin-transform of the hadronic cross section in xt = 2px/VS, with v^S" the 
hadronic cm. energy. The threshold logarithms exponentiate in transform space. 

Regarding phenomenology, the larger xt, the more dominant will the threshold logarithms be, 
and hence the more important will threshold resummation effects be. In addition, because of the 
convoluted form of the partonic cross sections and the parton distribution functions (PDFs), the 
steep fall-off of the PDFs with momentum fraction x automatically enhances the contributions 
from the threshold regime to the cross section, because it makes it relatively unlikely that the 
initial partons have very high cm. energy. This explains why partonic threshold effects often 
dominate the hadronic cross section even at not so high xt- Studies of cross sections for [8] 
pp — > hX (with h a high-p^ hadron) and [91 [101 [HI [12] pp — * 7A in the fixed-target regime, 
where typically 0.2 < xr < 0.7, indeed demonstrate that threshold-resummation effects dominate 
there and can be very large and important for phenomenology. They enhance the theoretical cross 
section with respect to fixed-order calculations. 

These observations suggest to study the resummation also for jet production at hadron col- 
liders, in particular when xt is rather large. An application of particular interest is the jet cross 
section at very high transverse momenta (pr ~ several hundreds GeV) at the Tevatron [T3J |T%] . 
for which initially an excess of the experimental data over next-to-leading order (NLO) theory 
was reported, which was later mostly attributed to an insufficient knowledge of the gluon distri- 
bution [15] . Similarly large values of xt are now probed in pp collisions at RHIC, where currently 
yfs = 200 GeV and jet cross section measurements by the STAR collaboration are already ex- 
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tending to pr > 40 GeV [H] . In both these cases, one does expect threshold resummation effects 
to be smaller than in the case of related processes at similar xt in the fixed-target regime, just 
because (among other things) the strong coupling constant is smaller at these higher p?. On the 
other hand, as we shall see, the effects are still often non-negligible. 

Apart from addressing these interesting phenomenological applications, we believe we also 
improve in this paper the theoretical framework for threshold resummation in jet production. 
There has been earlier work in the literature on this topic [U [T7J IE] ■ In Ref . [3] the threshold 
resummation formalism for the closely related dijet production at large invariant mass of the 
jet pair was developed to next-to-leading logarithmic (NLL) order. In [17], these results were 
applied to the single-inclusive jet cross section at large transverse momentum, making use of 
earlier work [B] on the high-p^ prompt-photon cross section, which is kinematically similar. As 
was emphasized in [4], there is an important subtlety for the resummed jet cross section related to 
the treatment of the invariant mass of the jet. The structure of the large logarithmic corrections 
that are addressed by resummation depends on whether or not the jet is assumed to be massless 
at partonic threshold, even at the leading-logarithmic (LL) level. This is perhaps surprising at 
first sight, because one might expect the jet mass to be generally inessential since it is typically 
much smaller than the jet's transverse momentum px and in fact vanishes for lowest-order partonic 
scattering. However, the situation can be qualitatively understood as follows [4]: let us assume 
that we are defining the jet cross section from the total four-momentum deposited in a cone of 
aperture R Q. Considering for simplicity the next-to-leading order, we can have contributions by 
virtual 2 — 2 diagrams, or by 2 — ► 3 real-emission diagrams. For the former, a single particle 
produces the (massless) jet, in case of the latter, there are configurations where two particles in 
the final state jointly form the jet. 

Then, for a jet forced to be massless at partonic threshold, the contributions with two partons 
in the cone must either have one parton arbitrarily soft, or the two partons exactly collinear. 
The singularities associated with these configurations cancel against analogous ones in the virtual 
diagrams, but because the constraint on the real-emission diagrams is so restrictive, large double- 
and single-logarithmic contributions remain after the cancellation. This will happen regardless 
of the size R of the cone aperture, implying that the coefficients of the large logarithms will be 
independent of R. These final-state threshold logarithms arising from the observed jet suppress 
the cross section near threshold. Their structure is identical to that of the threshold logarithms 
generated by the recoiling "jet", because the latter is not observed and is indeed massless at 
partonic threshold. The combined final-state logarithms then act against the threshold logarithms 
associated with initial-state radiation which are positive and enhance the cross section pQ. 

If, on the other hand, the jet invariant mass is not constrained to vanish near threshold, far 
more final states contribute- in fact, there will be an integration over the jet mass to an upper 
limit proportional to the aperture of the jet cone. As the 2 — > 3 contributions are therefore much 
less restricted, the cancellations of infrared and collinear divergences between real and virtual 
diagrams leave behind only single logarithms [4], associated with soft, but not with collinear, 
emission. Compared to the previously discussed case, there is therefore no double-logarithmic 
suppression of the cross section by the observed jet, and one expects the calculated cross section 
to be larger. Also, the single- logarithmic terms will now depend on the jet cone size R. 

t Details of the jet definition do not matter for the present argument. 
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The resummation for both these cases, with the jet massless or massive at threshold, has 
been worked out in [3]. The study [IT] of the resummed single- inclusive high-p^ jet cross section 
assumed massless jets at threshold. From a phenomenological point of view, however, we see no 
reason for demanding the jet to become massless at the partonic threshold. The experimental jet 
cross sections will, at any given p?, contain jet events with a large variety of jet invariant masses. 
NLO calculations of single-inclusive jet cross sections indeed reflect this: they have the property 
that jets produced at partonic threshold are integrated over a range of jet masses. This becomes 
evident in the available semi-analytical NLO calculations [191 l2~0l I2H [22] • For these, the jet cross 
section is obtained by assuming that the jet cone is relatively narrow, in which case it is possible 
to treat the jet definition analytically, so that collinear and infrared final-state divergences may be 
canceled by hand. This approximation is referred to as the "small-cone approximation (SCA)". 
Section II. E in the recent calculation in [21] explicitly demonstrates for the SCA that the threshold 
double-logarithms associated with the observed final-state jet cancel, as described above. 

In light of this, we will study in this work the resummation in the more realistic case of jets 
that are massive at threshold. We will in fact make use of the NLO calculation in the SCA 
approximation of [21] to "match" our resummed cross sections to finite (next-to-leading) order. 
Knowledge of analytical NLO expressions allows one to extract certain hard-scattering coefficients 
that are finite at threshold and part of the full resummation formula. These coefficients will be 
presented and used in our paper for the first time. 

We emphasize that the use of the SCA in our work is not to be regarded as a limitation to 
the usefulness of our results. First, the SCA is known to be very accurate numerically even at 
relatively large jet cone sizes of R ~ 0.7 [20J E3] • In addition, one may use our results to obtain 
ratios of the resummed over the NLO cross sections. Such "if-factors" are then expected to be 
extremely good approximations for the effects of higher orders even when one goes away from the 
SCA and uses, for example, a full NLO Monte-Carlo integration code that allows to compute the 
jet cross section for larger cone aperture and for other jet definitions (see, for example, Ref. |24|). 
We will therefore in particular present i^-factors for the resummed jet cross section in this paper. 

The paper is organized as follows: in Sec. [2] we provide the basic formulas for the single- 
inclusive-jet cross section at fixed order in perturbation theory, and discuss the SCA and the role 
of the threshold region. Section [3] presents details of the threshold resummation for the inclusive- 
jet cross section and describes the matching to the analytical expressions for the NLO cross section 
in the SCA. In Sec. H] we give phenomenological results for the Tevatron and for RHIC. Finally, 
we summarize our results in Sec. The Appendix collects the formulas for the hard-scattering 
coefficients in the threshold-resummed cross section mentioned above. 



2 Next-to-leading order single-inclusive jet cross section 

Jets produced in high-energy hadronic scattering, Hi(Pi)H 2 (P2) — ► ]et(Pj)X, are typically defined 
in terms of a deposition of transverse energy or four-momentum in a cone of aperture R in pseudo- 
rapidity and azimuthal-angle space, with detailed algorithms specifying the jet kinematic variables 
in terms of those of the observed hadron momenta [TBI [25l [26| [27] [28l [29\. QCD factorization 
theorems allow to write the cross section for single-inclusive jet production in hadronic collisions in 
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terms of convolutions of parton distribution functions with partonic hard-scattering functions |3Uj : 
da = ^2 dxidx 2 f a /H 1 {xi,fi%) h/H 2 (x 2 ,fi F ) da ab (xiP 1 ,x 2 P2,Pj, Hf, Hr) , (1) 

a,b J 

where the sum runs over all initial partons, quarks, anti-quarks, and gluons, and where \ip and \Lr 
denote the factorization and renormalization scales, respectively. It is possible to use perturbation 
theory to describe the formation of a high-py jet, as long as the definition of the jet is infrared- 
safe. The jet is then constructed from a subset of the final-state partons in the short-distance 
reaction ah — > partons, and a "measurement function" in the da a b specifies the momentum Pj of 
the jet in terms of the momenta of the final-state partons, in accordance with the (experimental) 
jet definition. 

The computation of jet cross sections beyond the lowest order in perturbative QCD is rather 
complicated, due to the need for incorporating a jet definition and the ensuing complexity of the 
phase space, and due to the large number of infrared singularities of soft and collinear origin at 
intermediate stages of the calculation. Different methods have been introduced that allow the 
calculation to be performed largely numerically by Monte-Carlo "parton generators", with only 
the divergent terms treated in part analytically (see, for example, Ref. |24j). 

A major simplification occurs if one assumes that the jet cone is rather narrow, a limit known as 
the "small-cone approximation (SCA)" [T9], [201 [2U [22] . In this case, a semi-analytical computation 
of the NLO single-inclusive jet cross section can be performed, meaning that fully analytical 
expressions for the partonic hard-scattering functions da a b can be derived which only need to be 
integrated numerically against the parton distribution functions as shown in Eq. (Trj). The SCA 
may be viewed as an expansion of the partonic cross section for small S = RJ cosh?], where 77 is 
the jet's pseudo-rapidity. Technically, the parameter S is the half-aperture of a geometrical cone 
around the jet axis, when the four-momentum of the jet is defined as simply the sum of the four- 
momenta of all the partons inside the cone [191 [21]. At small 5, the behavior of the jet cross section 
is of the form A\og(S) + B + 0(5 2 ), with both A and B known from Refs. [T91[2I]. Jet codes based 
on the SCA have the virtue that they produce numerically stable results on much shorter time 
scales than Monte-Carlo codes. Moreover, as we shall see below, the relatively simple and explicit 
results for the NLO single-inclusive jet cross section obtained in the SCA are a great convenience 
for the implementation of threshold resummation, particularly for the matching needed to achieve 
full NLL accuracy. 

It turns out that the SCA is a very good approximation even for relatively large cone sizes of up 
to R ~ 0.7 [2T1 [201 123], the value used by both Tevatron collaborations. Figured] shows comparisons 
between the NLO cross sections for single-inclusive jet production obtained using a full Monte- 
Carlo code [21] and the SCA code of [21], for pp collisions at cm. energy y/S = 1800 GeV and very 
high pt- Throughout this paper we use the CTEQ6M [31] NLO parton distribution functions. We 
have chosen two different jet definitions in the Monte-Carlo calculation. One uses a conventional 
cone algorithm [25]; the other the CDF jet definition [13J. One can see that the differences with 
respect to the SCA are of the order of only a few per cent. We note that similar comparisons 
in the RHIC kinematic regime have been shown in [21]. In their recent paper |16| . the STAR 
collaboration used R = 0.4, for which the SCA is even more accurate. 

Encouraged by this good agreement, we will directly use the SCA analytical results when 
performing the threshold resummation. As stated in the Introduction, this is anyway not a 
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Figure 1: Ratio between NLO jet cross sections at Tevatron at y S — 1800 GeV, computed with 
a full Monte-Carlo code [21] and in the SCA. The solid line corresponds to the jet definition 
implemented by CDF [13] (with parameter R sep = 1.3), and the dashed one to the standard cone 
definition [25] ■ in both cases the size of the jet cone is set to R = 0.7, and the CTEQ6M NLO 
parton distributions, evaluated at the factorization scale fip — Pt, were used. 



limitation, because we will also always provide the ratio of resummed over NLO cross sections 
( K- factors) , which may then be used along with full NLO Monte-Carlo calculations to obtain 
resummed cross sections for any desired cone size or jet algorithm. 

A further simplification that we will make is to consider the cross section integrated over all 
pseudo-rapidities of the jet. As was discussed in [8], this considerably reduces the complexity 
of the resummed expressions. By simply rescaling the resummed prediction by an appropriate 
ratio of NLO cross sections one can nonetheless obtain a very good approximation also for the 
resummation effects on the non-integrated cross section, at central rapidities [llj. To perform the 
NLL threshold resummation for the full rapidity-dependence of the jet cross section remains an 
outstanding task for future work. 

From Eq. (pQ), we find for the single- inclusive jet cross section integrated over all jet pseudo- 
rapidity r], in the SCA: 



p% da SCA (x T ) 



^2 / dXl f a / H l ( Xl > Pf) / dx 2fb/H 2 {x2,V 2 F ) 

„ u Jo Jo 



dp T 

a,b 

dx T 5[x T - / d-q — —j— , (2) 

V \J%i%2 J Jfi_ 2 dx^dr] 

where as before xt = 2p T / \fS is the customary scaling variable, and xt = %Pt/ V$ with s = X\X 2 S 
is its partonic counterpart, r) is the partonic pseudo-rapidity, r) = rj — \ In (x l/a^), which has the 
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limits fj + = —f}_ = ln[(l + y 1 — Xj)/xt\. The dependence of the partonic cross sections on fip 
and fiR has been suppressed for simplicity. The perturbative expansion of the da a b in the coupling 
constant as(nn) reads 

da ab (x 2 T , fj, R) =al(fj, R ) daf b \x%, fj) + a s (fi R ) da^(x%, fj, R) + ©(ajj) . (3) 

As indicated, the leading-order (LO) term da ab has no dependence on the cone size R, because 
for this term a single parton produces the jet. The analytical expressions for the NLO terms da^f 
have been obtained in [T^l [21] . It is customary to express them in terms of a different set of 
variables, v and w, that are related to xt and fj by 

x 2 T = Avw(l-v) , e 2 ^ = ^^. (4) 

1 — v 

Schematically, the NLO corrections to the partonic cross section for each scattering channel then 
take the form 

sda„ h \w,v,R) . . _,. . . /ln(l — u>) 

A a6 (u, i2) 5(1 - w) + B ab (v, R) 1 



aw at> V 1 — w 



+C ab {v,R) ^— —J + F ab (w,v,R), (5) 

where the "plus" -distributions are defined as usual by 

1 dwf(w)[g(w)}+ = [ dw{f{w) - f(l))g(w) , (6) 
Jo 

and where the F ab (w, v, 5) collect all terms without distributions in w. Partonic threshold corre- 
sponds to the limit w — *■ 1. The "plus" -distribution terms in Eq. (jSJ) generate the large logarithmic 
corrections that are addressed by threshold resummation. At order k of perturbation theory, the 

leading contributions are proportional to ctg y J • Performing the integration of these 

terms over f), they turn into contributions oc ctg ln 2fc (1 — x\), as we anticipated in the Introduc- 
tion, and as we shall show below. Subleading terms are down by one or more powers of the 
logarithm. We will now turn to the NLL resummation of the threshold logarithms. 



3 Resummed cross section 



The resummation of the soft-gluon contributions is carried out in Mellin-iV moment space, where 
they exponentiate [H [21 El IU El El E] ■ ^ ^ ne same time, in moment space the convolutions between 
the parton distributions and the partonic subprocess cross sections turn into ordinary products. 
For our present calculation, the appropriate Mellin moments are in the scaling variable xt B 

a{N)= Cdx 2 T (4)(*-9 gg^M . (7) 



*We drop the superscript "SCA" from now on. 
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In Mellin-iV space the QCD factorization formula in Eq. ([2]) becomes 

a.b 



(8) 



where the f^, H are the moments of the parton distribution functions, 



fa/H{^F)= I dxX N l f a /H{x,l4) 
JO 



and where 



dw / dv [4v(l — v)w] 
u Jo 



N+1 sda ab (w,v) 
dw dv 



(9) 



(10) 



The threshold limit w — > 1 corresponds to iV — > oo, and the LL soft-gluon corrections contribute 
as a^ln 2 ™ iV. The large- N behavior of the NLO partonic cross sections can be easily obtained by 
using 



[4.(1 -.)f +1 f(v) 



[4.(1 - v)] N+1 ( / ( £ ) + O 



Up to corrections suppressed by 1/JV it is therefore possible to perform the .-integration of the 
partonic cross sections by simply evaluating them at . = 1/2. According to Eq. when . = 1/2 
is combined with the threshold limit w = 1, one has = 1- It is worth mentioning that in 
the same limit one has rj = 0, and therefore the coefficients for the soft-gluon resummation for 
the rapidity-integrated cross section agree with those for the cross section at vanishing partonic 
rapidity. This explains why generally the resummation for the rapidity-integrated hadronic cross 
section yields a good approximation to the resummation of the cross section integrated over only 
a finite rapidity interval, as long as a region around r\ = is contained in that interval |llj . 

The resummation of the large logarithms in the partonic cross sections is achieved by showing 
that they exponentiate in Sudakov form factors. The resummed cross section for each partonic 
subprocess is given by a formally rather simple expression in Mellin-iV space: 



a 



(res) 



(N) = J2Ca b A%A b N J«J d N 

c.d 



a (int)afi- 
+cd N 



•cd 



a 



(12) 



where the first sum runs over all possible final state partons c and d, and the second over all possible 
color configurations I of the hard scattering. Except for the Born cross sections <x^^ (which we 
have presented in earlier work [8]), each of the iV-dependent factors in Eq. (|12p is an exponential 
containing logarithms in N. The coefficients C a b collect all iV-independent contributions, which 
partly arise from hard virtual corrections and can be extracted from comparison to the analytical 
expressions for the full NLO corrections in the SCA. Finally, the G^^ cd are color weights obeying 
Yli ^ib^cd = 1- The color interferences expressed by the sum over / appear whenever the number 
of partons involved in the process at Born level is larger than three, as it is the case here[§. Figure [2] 
gives a simple graphical mnemonic of the structure of the resummation formula, and of the origin 
of its various factors, whose expressions we know present. 



In a more general case without rapidity integration, the terms G^ b ^ C(i x o^™"^(A0 should be replaced by the 

(ao mum. 



color-correlated Born cross sections a^^P 
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Figure 2: Pictorial representation of the resummation formula in Eq. ( fl2i) . 



Effects of soft-gluon radiation collinear to the initial-state partons are exponentiated in the 
functions A^f, which read: 



lnA^= / / -\A a (a s (q 2 )) , (13) 

Jo 1 ~ z J^ F T 

(and likewise for b) where Q 2 = 2p\-.. Jfj is the exponent associated with collinear, both soft and 
hard, radiation in the unobserved recoiling "jet", 



In 4 



1 z N ~ 1 - 1 r f {1 ~ z)Q dq 2 A . . 2vs 1 



l-z 



A d (a s (q 2 )) + -B d (a s ((l - z)Q 2 )) . (14) 

(l-z) 2 Q 2 Q 1 



,2 



The function J'fi describes radiation in the observed jet. As we reviewed in the Introduction, 
this function is sensitive to the assumption made about the jet's invariant mass at threshold [I] 
even at LL level. Our choice is to allow the jet to be massive at partonic threshold, which is 
consistent with the experimental definitions of jet cross sections and with the available NLO 
calculations. In this case, J'§ is given as 

r 1 /-i _ i 

In J'§ = / — C' c {a s {{l - z) 2 Q 2 )) (jet massive at threshold) . (15) 

Jo 1 — z 

A similar exponent was derived for this case in j3]. The expression given in Eq. ffTSI) agrees with 
the one of [1] to the required NLL accuracy. Notice that J'^ contains only single logarithms, which 
arise from soft emission, whereas logarithms of collinear origin are absent. This is explicitly seen 
in the NLO calculations in the SCA [201 EI] in which there is an integration over the jet mass up 
to a maximum value of 0(5 ~ R), even when the threshold limit is strictly reached. Collinear 
contributions that would usually generate large logarithms are actually "regularized" by the cone 
size 8 and give instead rise to log(<5) terms in the perturbative cross sections. 

If, however, the jet is forced to be massless at partonic threshold the jet-function is identical 
to the function for an "unobserved" jet given in Eq. (fl4l) [U ITT] : 

In J'£ = In J^r (jet massless at threshold) , (16) 
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which produces a (negative) double-logarithm per emitted gluon, because it also receives collinear 
contributions due to the stronger restriction on the gluon phase space. There is then no dependence 
on log(<5) in this case. 

Finally, large-angle soft-gluon emission is accounted for by the factor A^' )a6 ^ c<i , which reads 

InA'^^ Z — -Z^ crf (a s ((l-^) 2 Q 2 )), (17) 

Jo i- — z 

and depends on the color configuration I of the participating partons. 

The coefficients A a , B a , C' a , Dja^d in Eqs. (TT3l) . (|T4l) . (fT5|) . (jT7|) are free of large logarithmic 
contributions and are given as perturbative series in the coupling constant 0:3: 



^(a s ) = ^ (1) +(-) 2 ^ (2) + --- 

71 V 71 / 



for each of them. For the resummation to NLL accuracy we need the coefficients Aa , A a 2 \ 
B a , C'a 1 , and D^ b ^ cd . The last of these depends on the specifics of the partonic process under 
consideration; all the others are universal in the sense that they only distinguish whether the 
parton they are associated with is a quark or a gluon. The LL and NLL coefficients A a , A^ and 
B a 1 ^ are well known 



A® = C a , A^ = \cjC, B^= la (19) 

with 

where C g = C A = N c = 3, C q = C F = (iV c 2 - l)/2N c = 4/3, 7, = -3/2C F = -2, lg = -2nb , and 
Nf is the number of flavors. The coefficients C' a needed in the case of jets that are massive at 
threshold, may be obtained by comparing the first-order expansion of the resummed formula to 
the analytic NLO results in the SCA. They are also universal and read 

CW = -C a ]og(?) . (21) 



This coefficient contains the anticipated dependence on log(5) that regularizes the final-state 
collinear configurations. As expected from a term of collinear origin, the exponent ffT5]) hence 
provides one power of log(<5) for each perturbative order. 

The coefficients Df^^^ governing the exponentiation of large-angle soft-gluon emission to NLL 
accuracy, and the corresponding "color weights" Gi a b^ c d, depend both on the partonic process 
and on its "color configuration" . They can be obtained from [H [T71 E31 El] where the soft 
anomalous dimension matrices were computed for all partonic processes. The results are given 
for the general case of arbitrary partonic rapidity 7). As discussed above, the coefficients for the 
rapidity- integrated cross section may be obtained by setting 17 = 0. We have presented the full 
set of the Dfl b ^ cd and Gi a b^ c d m the Appendix of our previous paper [8]. 

Before we continue, we mention that one expects that a jet cross section defined by a cone 
algorithm will also have so-called "non-global" threshold logarithms [33| [55] . These logarithms 
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arise when the observable is sensitive to radiation in only a limited part of phase space, as is the 
case in presence of a jet cone. For instance, a soft gluon radiated at an angle outside the jet cone 
may itself emit a secondary gluon at large angle that happens to fall inside the jet cone, thereby 
becoming part of the jet [35l [36]. Such configurations appear first at the next-to-next-to-leading 
order, but may produce threshold logarithms at the NLL level. One may therefore wonder if our 
NLL resummation formulas given above are complete. Fortunately, as an explicit example in [35] 
shows, it turns out that these effects are suppressed as R\og(R) in the SCA. They may therefore 
be neglected at the level of approximation we are making here. Given their only mild suppression 
as R — > 0, a study of non-global logarithms in hadronic jet cross sections is an interesting topic 
for future work. 



Returning to our resummed formulas, it is instructive to consider the structure of the leading 
logarithms. The LL expressions for the radiative factors are 



a a 

N 



J 



N 



exp 
exp 



L fOJn 2 (iV) 



[22) 



As discussed above, J' c does not contribute at the double-logarithmic level. Therefore, for a given 
partonic channel, the leading logarithms are 



a 



X £LW « exp 



(23) 



The exponent is positive for each partonic channel, implying that the soft-gluon effects will increase 
the cross section. This enhancement arises from the initial-state radiation represented by the A 0j & 
and is related to the fact that finite partonic cross sections are obtained after collinear (mass) 
factorization [TJ [2]. In the MS scheme such an enhancing contribution is (for a given parton 
species) twice as large as the suppressing one associated with final-state radiation in J d , for which 
no mass factorization is needed. For quark or anti-quark initiated processes, the color factor 
combination appearing in Eq. (I23p ranges from 2Cp — Cp/2 = 2 for the qq — > qq channel to 
2Cp — Ca/2 = 7/6 for qq — > gg, while for those involving a quark-gluon initial state one has larger 
factors, C F + C A - C F /2 = 11/3 (for qg -> qg) or C F + C A - C A /2 = 17/6 (for qg -> gq). Yet 
larger factors are obtained for gluon-gluon scattering, with 2C A — C A /2 = 9/2 for gg — ► gg and 
2C A — Cf/2 = 16/3 for gg — > qq. Initial states with more gluons therefore are expected to receive 
the larger resummation effects. We mention that if the observed jet is assumed strictly massless 
at threshold, an extra suppression term proportional to J c arises (see Eq. f [TB]) ). 

It is customary to give the NLL expansions of the Sudakov exponents in the following way [2]: 

lnA^asG4),Q7/4Q7/4) = InNh^(X) + / i ( 2 )(A,Q 2 //4;Q 2 / / 4) + {a s (a s lnN) k ) , (24) 
lnJ^(a s G4),Q 2 //4) = lniV/i 1 )(A) + /i 2 )(A,g 2 /^) + 0(« s (a s lniV) fc ) , (25) 



lnj;> s (/4)) 



In A 



(int) ab— >cd 
IN ' 



dp 

2nb 



<*s(/4)) 



ln(l-2A) + 0(as(a s lnAO fc ) , 

o 

lab— *cd o\N i /„ I „ l^AAfc" 1 



2tt6 



ln(l-2A) + 0(a s (a s \nN) h ) 



(26) 
(27) 
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with A = b as(fi 2 R ) \nN. The LL and NLL auxiliary functions hi 1 ' 2 ^ and fa 1 ^' are 



£•(1,2) 



^ 1} (A)= + 



+ 



2tt6 A 

4(2) 
2^2 

2nb 3 



[2A 



2A)ln(l -2A)] 



(2f 



4(1) 

[2A + ln(l-2A)]-^r^ln(l 



nb 



2A) 



1 



2A + ln(l-2A) + -ln 2 (l-2A) 
2 



4> g 2 

- T" A In — , 

7ib [Xp 



/i 1} (A) 
/i 2) (A,Q 2 //4) 



Aa 



(1) 



2tt6oA 
27r&jj 



A {1) (l 2 
+ ^-[2A + ln(l-2A)]ln^ 

27T&0 /X^ 



[1 - 2A) ln(l - 2A) - 2(1 - A) ln(l - A) 



ln(l - 2A) - 2 ln(l - A) + - ln 2 (l - 2A) - ln 2 (l - A) 



(29) 
(30) 



+ 



2tt6o 
A? 
2n 2 b 2 



ln(l - A) 



4% 



wb 

21n(l - A) - ln(l - 2A) 



ln(l - A) - ln(l - 2A) 

21n(l - A) -ln(l-2A) 



(31) 



AP 
2nb 



1 

In — , 

/4 



where bo, b\ are the first two coefficients of the QCD /^-function: 



1 



12tt 



(UC A -2N f ) 



1 



24tp 



(17C\ - 5C A N f - 3C F N f ) . 



(32) 



The ^-independent coefficients C a b in Eq. (|T2|) . which include the hard virtual corrections, 
have the perturbative expansion 



a 



„ /( = l + ^C7«+0(^). 

7T 



(33) 



The Cffc we need to NLL are obtained by comparing the (9(as)-expansion (not counting the 
overall factor a| of the Born cross sections) of the resummed expression with the fixed-order 
NLO result for the process, as given in [T9| [21]. The full analytic expressions for the C^J are 
rather lengthy and will not be given here. For convenience, we present them in numerical form 
in the Appendix. We note that apart from being useful for extracting the coefficients Ca 1 ^ and 
, the comparison of the (9(as)-expanded resummed result with the full NLO cross section also 
provides an excellent check of the resummation formula, since one can verify that all leading and 
next-to-leading logarithms are properly accounted for by Eq. (fl2"l) . 

The improved resummed hadronic cross section is finally obtained by performing an inverse 
Mellin transformation, and by properly matching to the NLO cross section da^ NLO \xT) / dpr as 
follows: 



p% da imatch) (x T ) 
dp T 



E 



Cmp+joo 



dN 



C M p-i°° 



2m 



: (4r N+1 f a %M) h%M) 



X 



a 



(res) 
ab^cd 



(AO 



a 



(res) 
ab— *cd 



(AO 



(NLO) 



+ 



plda^°\x T ) 
dp T 



(34) 
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where <3"^!!, C(i is given in Eq. (jl2~!) and (<r^!!; cd )(NLO) represents its perturbative truncation at NLO. 
Thus, as a result of this matching procedure, in the final cross section in Eq. ( 1341) the NLO cross 
section is exactly taken into account, and NLL soft-gluon effects are resummed beyond those 
already contained in the NLO cross section. 

The functions h a ' (A) and f a ' (A) in Eqs. (T28T )- fT3TT) are singular at the points A = 1/2 and/or 
A = 1. These singularities are related to the divergent behavior of the perturbative running 
coupling as near the Landau pole, and we deal with them by using the Minimal Prescription 
introduced in Ref. [2]. In the evaluation of the inverse Mellin transformation in Eq. flMl) . the 
constant Cmp is chosen in such a way that all singularities in the integrand are to the left of the 
integration contour, except for the Landau singularities, that are taken to lie to its far right. We 
note that an alternative to such a definition one could choose to expand the resummed formula 
to a finite order, say, next-to-next-to- leading order (NNLO), and neglect all terms of yet higher 
order. This approach was adopted in Ref. [T7j. We prefer to keep the full resummed formula in 
our phenomenological applications since, depending on kinematics, high orders in perturbation 
theory may still be very relevant [8] . It was actually shown in [2] that the results obtained within 
the Minimal Prescription converge asymptotically to the perturbative series. 

This completes the presentation of all ingredients to the NLL threshold resummation of the 
hadronic single-inclusive jet cross section. We will now turn to some phenomenological applica- 
tions. 



4 Phenomenological Results 

We will study the effects of threshold resummation on the single-inclusive jet cross section in pp 
collisions at = 1.8 TeV and y/S = 630 GeV cm. energies, and in pp collisions at v^S* = 
200 GeV. These choices are relevant for comparisons to Tevatron and RHIC data, respectively. 
Unless otherwise stated, we always set the factorization and renormalization scales to fiF = = 
Pt and use the NLO CTEQ6M [31] set of parton distributions, along with the two-loop expression 
for the strong coupling constant cts- 

We will first analyze the relevance of the different subprocesses contributing to single-jet pro- 
duction. The left part of Fig. [3] shows the relative contributions by "gg" (gg, qq', qq and qq' 
combined), qg and gg initial states at Born level (dashed lines) and for the NLL resummed case 
(without matching, solid lines). Here we have chosen the case of pp collisions at v^S* = 1.8 TeV. 
As can be seen, the overall change in the curves when going from Born level to the resummed case 
is moderate. The main noticeable effect is an increase of the relative importance of processes with 
gluon initial states toward higher p?, compensated by a similar decrease in that of the qq channels. 
In the right part of Fig. [3] we show the enhancements from threshold resummation for each initial 
partonic state individually, and also for their sum. At the higher pt, where threshold resummation 
is expected to be best applicable, the enhancements are biggest for the gg channel, followed by 
the qg one. All patterns observed in Fig. [3] are straightforwardly understood from Eq. (1231 . which 
demonstrates that resummation yields bigger enhancements when the number of gluons in the 
initial state is larger. We note that results very similar to those shown in the figure hold also at 
y/S = 630 GeV, if the same value of xt = 2p T /\/S is considered. This remains qualitatively true 
even when we go to pp collisions at y/S = 200 GeV, except for the larger enhancement in the 
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1.0 rr-T- 




100 200 300 400 100 200 300 400 



p T (GeV) P T (GeV) 

Figure 3: Left: relative contributions of the various partonic initial states to the single-inclusive jet 
cross section in pp collisions at V~S = 1.8 TeV, at Born level (dashed) and for the NLL resummed 
case (solid). We have chosen the jet cone size R = 0.7. Right: ratios between resummed and Born 
contributions for the various channels, and for the full jet cross section. 




quark contribution, due to the dominance of the qq channel (instead of the qq as in pp collisions) 
with a larger color factor combination in the Sudakov exponent (see the discussion after Eq. (T22]) V 

Figure H] repeats the studies made for Fig. [3] for this case. As one can see, if the same xt 
as in Fig. [3] is considered, the qq scattering contributions are overall slightly less important. At 
the same time, resummation effects are overall somewhat larger because the pt values are now 
much smaller than in Fig. [3], so that the strong coupling constant that appears in the resummed 
exponents is larger. 
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Figure 5: Ratio between the expansion to NLO of the (unmatched) resummed cross section and 
the full NLO one (in the SCA), for pp collisions at t/s = 1.8 TeV (solid) and y/s = 630 GeV 
(dots), and for pp collisions at y/s = 200 GeV (dashed). 



Before presenting the results for the matched NLL resummed jet cross section and X-factors, 
we would like to identify the kinematical region where the logarithmic terms constitute the bulk 
of the perturbative corrections and subleading contributions are unimportant. Only in these is 
the resummation expected to provide an accurate picture of the higher-order terms. We can 
determine this region by comparing the resummed formula, expanded to NLO, to the full fixed- 
order (NLO) perturbative result, that is, by comparing the last two terms in Eq. (1341 . Figure [5] 
shows this comparison for both Tevatron energies and for the RHIC function of the 

"scaling" variable xt- As can be observed, the expansion correctly reproduces the NLO result 
within at most a few per cent over a region corresponding to px > 200 GeV for the higher Tevatron 
energy, and to pr > 30 GeV at RHIC. This demonstrates that, at this order, the perturbative 
corrections are strongly dominated by the terms of soft and/or collinear origin that are addressed 
by resummation. The accuracy of the expansion improves toward the larger values of the jet 
transverse momentum, were one approaches the threshold limit more closely. 

Having established the importance of the threshold corrections in a kinematic regime of interest 
for phenomenology, we show in Fig. [H]the impact of the resummation on the predicted single-jet 
cross section at y/S =1.8 TeV. NLO and NLL resummed results are presented, computed at three 
different values of the factorization and renormalization scales, defined by fip = fiR — (pt (with 
( = 1,2, 1/2). The most noticeable effect is a remarkable reduction in the scale dependence of 
the cross section. This observation was also made in the previous study [TTj . If, as customary, 
one defines a theoretical scale "uncertainty" by A = (<j(C = 0.5) — cr(C = 2))/a(C = 1), the 
improvement is considerable. While A lies between 20 and 25% at NLO, it never exceeds 8% for 
the matched NLL result. The inset plot shows the NLL f^-factor, defined as 

(res) _ da^/dp T 
K ~ da^°)/dpT ' 1 ' 

at each of the scales. The corrections from resummation on top of NLO are typically very mod- 
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erate, at the order of a few per cent, depending on the set of scales chosen. The higher-order 
corrections increase for larger values of the jet transverse momentum. These findings are again 
consistent with those of [17], even though more detailed comparisons reveal some quantitative 
differences that must be related to either the different choice of the resummed final-state jet func- 
tion in jTTj (see discussion in Sec. [3]), or to the fact that [UJ uses only a NNLO expansion of the 
resummed cross section. The main features of our results remain unchanged when we go to the 
Tevatron-run II energy of v^S* = 1.96 TeV, at which measured jet cross sections are now avail- 
able [37] • Quantitatively very similar results are also found for the lower Tevatron center-of-mass 
energy, as seen in Fig. [7J In the case of pp collisions at y/s = 200 GeV, presented in Fig. [8] a 
similar pattern emerges, even though the resummation effects tend to be overall somewhat more 
substantial here. 
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Figure 6: NLO and NLL results for the single-inclusive jet cross section in pp collisions at v^S" = 1.8 
TeV, for different values of the renormalization and factorization scales. We have chosen R = 0.7. 
The inset plot shows the corresponding ^-factors as defined in Eq. (15B1) . 



In Fig. [9] we analyze how the resummation effects build up order by order in perturbation 
theory. We expand the matched resummed formula beyond NLO and define the "partial" soft- 
gluon i^-factors as 

da( ^ )/dpT , 

which for n — 2, 3, . . . give the additional enhancement over full NLO due to the 0(a 2 s +n ) terms 
in the resummed formula^. Formally, K 1 = 1 and K°° = of Eq. ([35]). The results for 

.ft-2,3,4,00 are given in the figure, for the case oi pp collisions at v^S* = 1.8 TeV. One can see that 



^We recall that the Born cross sections are of 0(aa), hence the additional power of two in this definition. 
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Figure 7: Same as Fig. [6l but for v^S* =630 GeV. 

contributions beyond N 3 LO (n = 3) are very small, and that the 0(a^) result can hardly be 
distinguished from the full NLL one. 

It is interesting to contrast the rather modest enhancement of the jet cross section by resum- 
mation to the dramatic resummation effects that we observed in j8] in the case of single-inclusive 
pion production, H\H-2, — > 7rX, in fixed-target scattering at typical cm. energies of \fS ~ 30 GeV. 
Even though in both cases the same partonic processes are involved at Born level, there are several 
important differences. First of all, the values of pt are much smaller in fixed-target scattering 
(even though roughly similar values of xt = 2p T /\fS are probed), so that the strong coupling 
constant «s(pr) is larger and resummation effects are bound to be more significant. Furthermore, 
for the process H1H2 — > ttX one needs to introduce fragmentation functions into the theoretical 
calculation that describe the formation of the observed hadron from a final-state parton. As the 
hadron takes only a certain fraction z > 0.5 of the parent parton's momentum, the partonic hard- 
scattering necessarily has to be at the higher transverse momentum pr/z in order to produce a 
hadron with p^- Thus one is closer to partonic threshold than in case of a jet produced with 
transverse momentum px which takes all of a final-state parton's momentum. In addition, it 
turns out [H EH] that due to the factorization of final-state collinear singularities associated with 
the fragmentation functions, the "jet" function J'^ in the resummation formula Eq. (|T2j) is to be 
replaced by a factor A C N , which has enhancing double logarithms. 

Finally, as one illustrative example, we compare our resummed jet cross section to data from 
CDF [13J at \/S = 1800 GeV. While so far we have always considered the cross section integrated 
over all jet rapidities, we here need to account for the fact that the data cover only a finite 
region in rapidity, 0.1 < \r]\ < 0.7. Also, we would like to properly match the jet algorithm 
chosen in experiment, rather than using the SCA. We have mentioned before that both these 
issues can be accurately addressed by "rescaling" the resummed cross section by an appropriate 
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Figure 8: Same as Fig. El but for pp collisions at =200 GeV and R = 0.4. 



ratio of NLO cross sections. We simply multiply our i^-factors defined in Eq. (1551) and shown 
in Fig. [6] by da^ MC \0.1 < \t)\ < 0.7) /dpr, the NLO cross section obtained with a full Monte- 
Carlo code [23], in the experimentally accessed rapidity regime. The comparison between the 
data and the NLO and NLL-resummed cross sections is shown in Fig. [10] in terms of the ratios 
"(data— theory) /theory". As expected from Fig. 0, the impact of resummation is moderate and 
in fact smaller than the current uncertainties of the parton distributions [31] • Nonetheless, it does 
lead to a slight improvement of the comparison and, in particular, the plot again demonstrates 
the reduction of scale dependence by resummation. 



5 Conclusions 

We have studied in this paper the resummation of large logarithmic threshold corrections to 
the partonic cross sections contributing to single-inclusive jet production at hadron colliders. Our 
study differs from previous work [17] mostly in that we allow the jet to have a finite invariant mass 
at partonic threshold, which is consistent with the experimental definitions of jet cross sections 
and with the available NLO calculations. Moreover, using semi-analytical expressions for the 
NLO partonic cross sections derived in the SCA [T91 [2T], we have extracted the N— independent 
coefficients that appear in the resummation formula, and properly matched our resummed cross 
section to the NLO one. We hope that with these improvements yet more realistic estimates of 
the higher-order corrections to jet cross sections in the threshold regime emerge. 

It is well known that the NLO description of jet production at hadron colliders is overall very 
successful, within the uncertainties of the theoretical framework and the experimental data. From 
that perspective, it is gratifying to see that the effects of NLL resummation are relatively moderate. 
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On the other hand, resummation leads to a significant decrease of the scale dependence, and we 
expect that knowledge of the resummation effects should be useful in comparisons with future, 
yet more precise, data, and for extracting parton distribution functions. Given the general success 
of the NLO description, we have mostly focused on .fT-factors for the resummed cross section 
over the NLO one, and only given one example of a more detailed comparison with experimental 
data. We believe that these .ff -factors may be readily used in conjunction with other, more flexible 
NLO Monte-Carlo programs for jet production, to estimate threshold-resummation effects on cross 
sections for other jet algorithms and possibly for larger cone sizes. 
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Appendix: First-order coefficients C ah in the SCA 

In this appendix we collect the process-dependent coefficients Cj^ for the various partonic channels 
in jet hadroproduction in the SCA. The C$ are constant, that is, they do not depend on the 
Mellin-moment variable N. They may be extracted by expanding the resummed cross section in 
Eq. (j!2p to first order and comparing it to the full NLO cross section in the SCA. For the sake 
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Figure 10: Ratios "(data— theory) /theory" for data from CDF at 1.8 TeV [13] and for NLO and 
NLL resummed theoretical calculations. We have chosen the theory result at scale as "default" ; 
results for other scales are also displayed in terms of their relative shifts with respect to the default 
theory. 



of simplicity, we provide the only in numerical form, as the full analytic coefficients have 
relatively lengthy expressions. We find: 



Ml) 

^qq> 
qq' 

Mi) 

^qq 



Ml) 
^qq 

,7(1) 
M 1 ) 

^99 



17.9012 + 1.38763 log 
19.0395 + 1.38763 log 



R 
R 



R 

13.4171 + 1.6989 log-, 

R 

17.1973 + 1.38763 log - 

R 

14.4483 + 2.58824 log — 

R 

14.5629 + 3.67884 log - 



(37) 



where R is the jet cone size. 
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